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Ionized Gas in the First 10 kpc of the Interstellar Galactic Halo 

J. Christopher Howk^’^, Kenneth R. Sembach^, & Blair D. Savage"^ 

ABSTRACT 

We present Far Ultraviolet Spectroscopic Explorer (FUSE) observations of the post-asymptotic 
giant branch star von Zeipel 1128 {I = 42?5, b = -|-78?7; d = 10.2 kpc, z = 10.0 kpc; ulsr = 
—140± 8 km s“^), located in the globular cluster Messier 3. The FUSE observations cover the 
wavelength range 905 - 1187 A at ~ 20 km s“^ (FWHM) resolution. These data exhibit many 
photospheric and interstellar absorption lines, including absorption from ions associated with 
the warm neutral, warm ionized, and highly-ionized phases of the interstellar medium along this 
sight line. We derive interstellar column densities of H I, P II, Ar I, Fe II, Fe III, S III, and O VI, 
with lower limits for C II, C III, N I, O I, and Si II. Though the individual velocity components 
within the absorption profiles are unresolved by FUSE, a comparison of the velocity distribution 
of depleted or ionized species with the neutral species suggests the thick disk material along this 
sight line is infalling onto the Galactic plane, while material near the plane is seen closer to rest 
velocities. Ionized hydrogen represents > 12%, most likely ^ 45%, of the total hydrogen column 
along this sight line, most of it associated with the warm ionized phase. The warm ionized and 
neutral media toward von Zeipel 1128 have very similar gas-phase abundances and kinematics: 
the neutral and ionized gases in this region of the thick disk are closely related. Strong O VI 
absorption is seen with the same central velocity as the warm ionized gas, though the O VI 
velocity dispersion is much higher (cr = ^/2b = 32 km s“^). Virtually all of the O VI is found 
at velocities where lower-ionization gas is seen, suggesting the O VI and WNM/WIM probes 
are tracing different portions of the same structures (e.g., the O VI may reside in interfaces 
surrounding the WNM/WIM clouds). We see no evidence for interstellar absorption associated 
with the globular cluster Messier 3 itself nor with the circumstellar environment of von Zeipel 
1128. Neither high velocity cloud absorption (with |wlsr| ^ 125 km s“^) nor high velocity- 
dispersion gas (with tr ~ 60 km s“^) is seen toward von Zeipel 1128. 

Subject headings: ISM: atoms - ISM: structure - ultraviolet: ISM 


1. Introduction 

There is now substantial evidence that the interstellar medium (ISM) of the Milky Way is best described 
as a “multiphase” medium. In standard equilibrium models of the ISM (e.g., McKee & Ostriker 1977, Wolfire 
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et al. 1995), the term “phase” refers to gaseous regions of distinct temperatures, densities, and ionization 
states. Most models assume these regions are physically independent and exist in (sometimes rough) pressure 
equilibrium. Extraplanar material in spiral galaxies seems to exhibit all of the principal phases found in 
the disk of the Milky Way: the cold neutral medium [CNM; T ~ 10^ K; a;(H II) = n(H II)/n(Htot) ~ 10“^; 
Howk & Savage 1999a, 2000; Richter 2001]; the warm neutral medium [WNM; T ~ 10^ K; x(H II) ~ 0.1; 
Dickey & Lockman 1990]; the warm ionized medium ]WIM; a;(H II) > 0.8; Reynolds 1993]; and the hot ISM 
[T ~ 10® K; xin II) ~ 1.0; Snowden et al. 1998]. 

While the evidence strongly suggests their presence in the extraplanar regions of the Milky Way, the 
properties and physics of these phases and their relationships are as yet poorly characterized for extraplanar 
gas. As such, there are many unanswered questions regarding the nature of gas in the thick disk^ of the 
Milky Way and other spiral galaxies. What fraction of the interstellar thick disk is associated with each of 
these phases? What are the physical conditions within each of these phases? How does gas far from the 
midplane differ from that in the interstellar thin disk? What implications can be drawn from a comparison 
of the conditions of gas associated with the thin and thick disks? Are the identified phases simply artifacts 
of the manner in which the gas of the Galactic thick disk is studied? Similarly, are the various tracers of 
thick disk matter probing different portions of the same structures? 

While there is still relatively little information available to answer such questions, the relationship 
between the warm neutral and warm ionized gas of the thick disk is the best studied. Spitzer & Fitzpatrick’s 
(1993) study of the sight line to the distant star HD 93521 suggests a strong relationship between the warm 
neutral and ionized phases of the Galactic thick disk. They used high-resolution absorption line spectroscopy 
to show that excited-state C II (i.e., C 11*), produced by collisions of ground-state C II with warm electrons, 
exhibits the same velocity structure as other low-ionization species thought to trace the WNM (e.g., Fe II, 
Si II, and S II). Spitzer & Fitzpatrick argued that the electrons and the neutral hydrogen were well-mixed 
in a partially-ionized gas. Howk & Savage (1999b) found a similar kinematic alignment for the sight line 
toward p Leo (HD 91316), which lies near HD 93521 on the sky. Along this sight line absorption from 
low-ionization lines such as Zn II were shown to have a velocity structure very similar to that seen in the 
moderately-ionized species S III and Al III. The latter ions have ionization potentials significantly higher 
than that of H I, making them unlikely to occur in regions dominated by neutral hydrogen. 

Reynolds et al. (1995) approached this issue by comparing the distribution of Ha and H I 21-cm 
emission over a 10° x 12° region of the sky centered on {l,b) = (144°,—21°). They noted the existence 
of “Ha-emitting H I clouds,” structures seen in both the Ha and H I maps that share morphological and 
kinematical properties when viewed at the low spatial resolution of their typical maps. They showed the 
ratio of Ha intensity to H I column density increases with height above the plane in this direction. While 
Reynolds et al. found several prominent examples of structures that seem to bear both neutral and ionized 
hydrogen, when they examined one such filament at higher spatial resolution, they found that the neutral 
and ionized gases tended to be displaced from one another on the sky. While the Ha and H I were associated 
with the same filament, the ionized and neutral gases traced different portions of this structure. 

If these examples are representative of the high-latitude sky, they suggest much of the gas associated 


®In this paper we will often use the term “thick disk” to refer to the extraplanar distribution of interstellar matter in the 
Milky Way. This is analogous to the term “halo” used by many authors. We make this choice of terminology both because the 
ISM within the first few kiloparsecs of the midplane tends to rotate as a disk and because we wish to distinguish this lower-lying 
material from the much more extended corona, whose presence is inferred by observations of O VI associated with very distant 
gas clouds (e.g., high-velocity clouds; see Sembach et al. 2000, 2003). 
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with the neutral and ionized phases of the Galactic thick disk are distributed in a similar manner with 
similar physical conditions. If the Spitzer & Fitzpatrick proposition that neutral and ionized hydrogen are 
intermixed in the same structures is true, then the ionization of hydrogen (and the S III and A1 III toward 
p Leo) within these “clouds” most likely occurs by means other than photoionization by OB stars (e.g., 
shocks, X-ray photoionization, or cosmic-ray ionization), although it is generally thought that only ionizing 
radiation from OB stars has enough power to provide for the observed ionization of the WIM. Furthermore, 
the ionization of hydrogen through mechanisms such as shocks and X-ray photoionization may also produce 
by-products such as the highly-ionized gas traced by O VI and C IV. The work of Reynolds et al. (1995) 
suggests a close relationship between the warm neutral and ionized gases, though one in which these phases 
are indeed distinct. 

These studies are all difficult to generalize to the whole high-latitude sky. Too few sight lines have 
been studied in the detailed manner that the three discussed above have been to draw conclusions about the 
general relationship between the phases of the Galactic thick disk. Furthermore, these studies provide little 
insight into the role of hotter, more highly-ionized material. 

This paper presents Far Ultraviolet Spectroscopic Explorer (FUSE) observations of the post-asymptotic 
giant branch (post-AGB) star von Zeipel 1128 (hereafter vZ 1128; von Zeipel 1908). This star lies in the 
globular cluster Messier 3 (NGG 5272), located 10.2 kpc from the Sun at a height above the Galactic midplane 
of 10.0 kpc (Harris 1996). Thus, the sight line to vZ 1128 probes nearly all of the thick disk and much of 
the corona of the Galaxy, including the most distant reaches typically only accessible via observations of 
extragalactic targets (e.g.. Savage et al. 2000; Savage et al. 2003). 

We use the present FUSE observations to study the content and abundances of the ionized gas, as well 
as its kinematics, along this extended Galactic sight line. We will emphasize the ionization characteristics 
of the sight line and the connection between the various “phases” of the ISM in this direction. 

In §2 we discuss the properties of the star and sight line. We summarize the FUSE observations and our 
reduction of the data in §3, while we summarize our analysis methodology in §4. We present the properties 
of the highly-ionized gas along the vZ 1128 sight line in §5 and of the warm ionized and neutral media along 
this path in §6. In §7 we detail the lack of high-velocity absorption along this sight line. In §8 we discuss 
these results, and we summarize our principle conclusions in §9. 


2. Stellar and Interstellar Sight Line Properties 

Table 1 summarizes the properties of the star vZ 1128 and the interstellar path to this object. First 
cataloged by von Zeipel (1908), vZ 1128 (also known as NGG 5272 ZNG 1, K 728, 11-57) is classified as a 
hot (Teff « 33, 000 K; Dixon et al. 1994) post-AGB star. Such objects are thought to be the remnants of 
AGB stars that have ceased nucleosynthetic burning. By this stage of evolution, stars have shed significant 
amounts of mass (i.e., their complete outer envelopes) and are on their way to becoming white dwarfs. 
Not all post-AGB stars pass through a planetary nebula phase. Lower mass objects leave the AGB before 
the thermal pulsing phase of AGB evolution. These low-mass stars evolve very slowly (by the standards 
of higher-mass post-AGB stars); by the time they are hot enough to provide significant ionizing flux, the 
material that would have become a planetary nebula has since dispersed (see review by Moehler 2001). 

Relatively few observations of vZ 1128 have been reported in the literature. Studies of the stellar 
properties include works by de Boer (1985) using the International Ultraviolet Explorer (lUE) and Dixon et 
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al. (1994) using the Hopkins Ultraviolet Telescope. The stellar properties arising from these studies place 
vZ 1128 on theoretical tracks of post-AGB evolution (e.g., Schonberner 1983) with a mass of ~ 0.55 Mq 
(D ixon et al. 1994), though model atmosphere fits give slightly smaller values. The radial velocity reported 
for this star in Table 1 is from the present work. The quoted error is dominated by the uncertainties in the 
absolute velocity calibration of FUSE (discussed below). 

The interstellar sight line to vZ 1128 traces a pathlength of ~ 10.2 kpc almost perpendicular to the 
Galactic plane (Djorgovski 1993). This extended path samples gas in the northern half of the Galactic thin 
disk, thick disk, and extended corona. Given the high latitude of this sight line, the path through the thin 
disk of the Milky Way is quite short, approximately the half-thickness of this layer. Much of the absorption 
observed along this sight line is therefore expected to be from material significantly above the Galactic plane. 
The interstellar sight line to vZ 1128 has previously been discussed by de Boer & Savage (1984), who note 
the presence of negative velocity absorption in a high-resolution (low signal to noise) lUE spectrum. They 
comment that this material may be associated with the infall of matter participating in a “galactic fountain.” 

Figure 1 shows emission line profiles for ionized and neutral hydrogen in the general direction of vZ 1128. 
The top panel shows the Ha spectrum towards vZ 1128 from the Wisconsin H-alpha Mapper (WHAM; 
Haffner et al. 2002; Reynolds et al. 1998b) and represents an average of the eight pointings centered within 
1?5 of vZ 1128. WHAM has a 1° beam and 12 km s“^ velocity resolution. The next two panels show H I 
21-cm spectra from the Leiden-Dwingeloo Survey (LDS; Hartmann & Burton 1997) with a 30' beam and 
from the NRAO 140-ft telescope (unpublished from the Danly et al. 1992 survey) with a 21' beam. The 
LDS spectrum is centered ~ 12' from vZ 1128 at {l,b) = (42?5, -|-78?5). Both H I spectra have a velocity 
resolution of ~ 1 km s“^ and have been cleaned of stray radiation. The weak positive velocity wing (centered 
at ■clsr ~ +30 km s“^) in the LDS spectrum may result from imperfect baseline or stray radiation correction. 

These spectra demonstrate the presence of neutral and ionized material over the velocity range —70 < 
■^LSR ^ +30 km s“^. The total H I column densities over this range measured by the LDS and NRAO 
spectra are A^(H I) = 1.0 x 10^° and 1.1 x 10^° cm“^, respectively. An analysis of the Lyman-a absorption in 
recent Space Telescope Imaging Spectrograph echelle-mode data suggest A^(H I) = (9.3 ± 0.7) x 10^® cm“^ 
(Howk, Sembach, & Savage 2003). The difference in these column densities is likely a result of H I structure 
on scales smaller than the beam size of the radio observations. 

The H I emission spectra show a strong peak near wlsr ~ ~4 km s“^, which is close to the expected 
location of gas in the Galactic thin disk or material in the halo that is corotating with the disk. The next 
most prominent component or blend of components is centered near wlsr ~ ~25 km s“^. These two blends 
have nearly equal column densities and account for 99% of the total H I column in the NRAO spectrum. 

There is little evidence for significant amounts of cold material along this sight line in the H I spectrum. 
While Gaussian decomposition of the H I line profiles is at best non-unique, we have found that reasonable 
fits yield 6 « 10 km s“^ for the strongest component centered at vlsr ~ ~4 km s“^, where b is the Doppler 
parameter [b = (2fcT/m)^/^ for pure thermal broadening; Spitzer 1978]. This corresponds to a kinetic 
temperature of T < 6,000 K. (The LDS spectrum yields a slightly smaller 6-value: 6 « 7 km s“^.) The gas 
at more negative velocities can be fitted by Gaussians with widths 9 + 6 < 15 km s“^ implying upper limits 
on the temperature between 5,000 K and 14,000 K. While a small amount of cool gas could be blended 
within the emission profile, most of the neutral material along this sight line seems to be associated with the 
WNM. 


The negative-velocity gas is seen more prominently in the Ha spectrum than the H I spectra, which 
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in the simplest interpretation suggests this infalling material has a greater degree of ionization.® The total 
intensity of the Ha line averaged over the pointings centered within 1?5 is 0.39 Rayleighs, which corresponds 
to an emission measure of « 0.68 (T/6000)® ® pc cm“®. The intensity of the nearest one degree field to 
vZ 1128 is /(Ha) = 0.17 ± 0.03, or EM « [0.30 ± 0.06] (T/6000)® ® pc cm“®. If the mean density of the 
ionized material is assumed to be Ue ~ 0.08 (see Reynolds 1991a), then the H II column density along the 
sight line to vZ 1128 (averaged over the 1° beam) is A^(H II) ~ 1.2 x 10^® cm“^ for T = 6, 000 K (~ 1.8 x 10^® 
cm“^ for T = 10, 000 K). In this case, the column density of ionized gas is ~ 13% (~ 19%) that of the neutral 
gas (see also §6.2). The average ratio of the column densities 7V(H II)/iV(H I) integrated vertically through 
the Galaxy is between ~ 20% — 30% (Gomez, Benjamin, & Gox 2001; Reynolds 1993; Dickey & Lockman 
1990). We note there is large degree of variation in the Ha intensities in the WHAM measurements of this 
direction, with the beam centered closest to vZ 1128 having the lowest intensity of those pointings centered 
within 1?5. The ratio of ionized to neutral gas can be as high as 29% (44% if T = 10,000 K) assuming an 
emission measure averaged over the pointings centered within 1?5 (see Figure 1). These estimates of A^(H II) 
should be considered to have large uncertainties since the assumption of a single (poorly known) average 
electron density is unlikely to be correct. Identifying a pulsar in M3 would be extremely helpful to this 
investigation since it would allow the direct determination of N(Il II) through its dispersion measure. 

The direct comparison of these spectra, which sample large regions of the sky (~ 0?35 — 2?0), with our 
absorption line measurements is complicated by the difference in beam size of the observations. However, 
the H I and Ha measurements provide information on the average properties of the neutral and ionized 
material in this general direction. In particular they show evidence for significant departures from the 
expected Galactic rotation, since gas along high-latitude sight lines participating in this rotation should be 
seen nearly at rest relative to the Local Standard of Rest (LSR; see §6.4). A comparison of the H I and Ha 
profiles also reveals the increased importance of ionized gas in the negative-velocity material compared with 
that near the LSR. 

Hot gas (T > 10® K) in the Galactic halo can be traced through its soft X-ray emission (Snowden et 
al. 1997, 1998). The ROSAT 1/A-keY count rate averaged over a 36' region toward vZ 1128 is 1178 ± 142 
counts s“^ arcmin“^ (Savage et al. 2003). For comparison, the average 1/4-keV count rate for latitudes 
b > 60° derived from the maps of Snowden et al. (1997) is 1034 ± 255 counts s“^ arcmin“^, where the 
standard deviation is quoted as the uncertainty. While the interpretation of soft X-ray count rates is strongly 
dependent upon the total column density of absorbing material along the sight line, it seems the count rate 
along the sight line toward vZ 1128 is comparable to the high-latitude average. 


3. Observations and Rednctions 

Table 2 gives a log of the FUSE observations used for this work. The planning and in-orbit performance 
of FUSE are discussed by Moos et al. (2000) and Sahnow et al. (2000). FUSE consists of four co-aligned 
telescopes and Rowland-circle spectrographs feeding two microchannel plate (MGP) detectors with helical 
double delay line anodes. Two of the telescope/spectrograph channels have SiC coatings providing reflectivity 
over the range ~ 905 — 1105 A, while the other two have AhLiF coatings for sensitivity in the ~ 1000 — 1187 
A range. We refer to these as SiC and LiF channels, respectively. 


®It should be noted that while the H I spectrum is weighted by the neutral hydrogen density, the Ha spectrum is weighted 
by the square of the electron density. Furthermore, the Ho: spectrum shown is an average over a field of radius 2°. 
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The light reflected from each mirror passes through a focal plane assembly containing three entrance 
apertures; all of the vZ 1128 data were taken through the 30" x 30" (LWRS) apertures. The four holographically- 
ruled spherical gratings project the spectra onto two detectors (each made up of two segments - segments A 
and B), with a LiF and SiC channel projected in parallel onto each detector. For the vZ 1128 observations, 
the data were taken in time-tag mode, so that each detected photon is recorded by its position and arrival 
time. 

The individual exposures for the three sets of observations (see Table 2) were merged into a single 
spectrum by concatenating the individual photon event lists. We employed the CALFUSE (v2.0.5) pipeline for 
transforming the photon event lists to a calibrated one-dimensional spectrum. CALFUSE was used to screen 
the photon lists for valid data with constraints imposed for Earth limb angle avoidance and passage through 
the South Atlantic Anomaly. Corrections for detector backgrounds, Doppler shifts caused by spacecraft 
orbital motions, geometrical distortions, and astigmatism were applied (Sahnow et al. 2000). The CALFUSE 
processing is described in more detail by Oegerle et al. (2000). 

The processed data have a nominal spectral resolution of < 20 km s“^ (FWHM), with a Icr relative 
wavelength dispersion solution accuracy of < ±6 km s“^. The zero point of the wavelength scale for individual 
FUSE observations is poorly determined. For this work we compare the H I 21-cm line observations of Danly 
et al. (1992) with the absorption profiles of Ar I A1048.220 and A1066.660 to correct the FUSELiFlA data to 
the local standard of rest.^ When used as a model optical depth distribution and converted to an absorption 
line profile, the H I emission data are a reasonable approximation of the Ar I profiles. The other channels 
were then shifted to match the LiFlA LSR velocities. We estimate that our approach to determining the 
absolute velocity zero point for the EUSE observations is accurate to ~ ±5 km s“^. 

Figure 2 shows the full EUSE spectrum of vZ 1128, with prominent interstellar features marked. The 
displayed data are taken from several different detector channels, with the majority from SiClB, LiFlA, and 
LiF2A. Most of the unmarked absorption lines in Figure 2 are narrow stellar features. We have measured 
the centroids of several (10) of these (Table 3) and determine an average velocity of r^LSR = —140 ± 8 
km s“^. Our determination compares favorably with the determination of Strom & Strom (1970) who 
And wlsr = —142± 15 km s“^. Our quoted uncertainty includes contributions from the uncertainties in the 
absolute velocity calibration and the FUSE relative wavelength scale. However, the standard deviation in our 
individual measurements, cr = ±2.4 km s“^, is significantly smaller than the ±6 km s“^ relative wavelength 
calibration uncertainties typically quoted for FUSE data. Due to the nature of the FUSE detectors (see 
Sahnow et al. 2000), certain regions of the detector (certain wavelengths) are better At by the average 
dispersion solution than others. Thus, the ±6 km s“^ uncertainty quoted is a worst-case value. We note 
that consistent measurements of stellar velocities in wavelength regions near the O VI 1031.926 A and Ar I 
1066.660 A transitions imply the relative velocities of these two interstellar absorption lines are very well 
calibrated (to within a few km s“^). 

Figure 2 demonstrates an important benefit of using hot post-AGB stars to probe the ISM: the stellar 
continuum is extremely simple. In large part this is because post-AGB stars generally lack the prominent 
stellar P-Gygni profiles caused by the line-driven winds of most early-type stars. Furthermore, while many 
photospheric absorption lines are present, they are quite narrow (with breadths smaller then the ISM lines 
for vZ 1128) and at velocities that generally place them far from interstellar lines of interest (at least for 


^We adopt the standard definition of the LSR, assuming a solar motion of +20 km in the direction (q,5)i900 = 
(18^,+30°) [{l,b) ~ (56°,+23°)]. This gives ulsr = ^helio + 10.9 km s“l in the direction of vZ 1128. For comparison, the 
Mihalas & Binney (1981) definition of the LSR requires a shift of the velocities used here by —1.16 km s“^. 
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vZ 1128). Indeed, the spectrum of vZ 1128 is comparable in simplicity to many AGN spectra but is at least 
twice as bright as 3C273 (e.g., Sembach et al. 2001). 


4. Analysis of Interstellar Absorption Lines 


Each interstellar feature was normalized by fitting a low-order Legendre polynomial to the adjacent 
stellar continuum. For vZ 1128, this fitting process was for the most part unambiguous due to the lack 
of interfering stellar features (the features that are present are easy to distinguish and avoid in the fitting 
process). The normalized absorption profiles for several important interstellar transitions are shown in Figure 
3. The empirically-determined signal-to-noise ratios of these data are between ~ 15 and 25 per 20 km s“^ 
resolution element. 

We measured equivalent widths of interstellar transitions following Sembach & Savage (1992), including 
their treatment of the uncertainties. Equivalent width measurements are given in Table 4 for numerous 
interstellar transitions present in the FUSE data. The uncertainties contain both statistical (photon and 
fixed-pattern noise) and continuum placement uncertainties. For most transitions we quote equivalent widths 
measured in detector 1 and 2 data separately. For wavelengths A > 1000 A the measurements are from LiFl 
and LiF2 data, for A < 1000 A from SiCl and SiC2 data. The velocity ranges over which the equivalent 
width integrations were carried out are given in the last column of Table 4. The adopted atomic data are 
also summarized in this table. The rest wavelengths and most of the /-values are drawn from the revised 
compilation of D. Morton (1999, private communication), except where noted below. 

Table 4 also gives the integrated apparent column densities (Savage & Sembach 1991) for each of the 
transitions using the atomic data given in the second and third columns. The apparent optical depth, Ta(v), 
is an instrumentally-blurred version of the true optical depth of an absorption line, given by 


ra(v) = - In [I(v)/Ic(v)] 


( 1 ) 


where Ic(v) is the estimated continuum intensity and I(v) is the observed intensity of the line as a function of 
velocity. This is related to the apparent column density per unit velocity, Na{v) [atoms cm“^ (km s“^)“^], 

by 


Na{v) 


meCTa{v) 

7re^ /A 


= 3.768 X lO^"* 



( 2 ) 


where A is the wavelength in A, and / is the atomic oscillator strength. In the absence of unresolved 
saturated structure the Na{v) profile of a line is a valid, instrumentally-blurred representation of the true 
column density distribution as a function of velocity, N(v). In the presence of unresolved saturated structure, 
the Na{v) distribution is always a lower limit to the true column density distribution. 


Table 5 gives our adopted interstellar column densities for the vZ 1128 sight line. Two methods are used 
to determine these values and limits. For cases where it is expected that the line profiles are fully resolved 
or that no saturation corrections need be made, we make use of the apparent column density measurements 
from Table 4 to determine the column densities. In the case of interstellar O VI absorption, the absorption 
profiles are typically broad enough to be fully resolved by the FUSE spectrographs. For the vZ 1128 sight 
line, this is demonstrated by the fact that the weak and strong members of the doublet give consistent 
apparent column densities (see Table 4). For other species, this is not necessarily the case. In cases where no 
saturation correction is required we adopt the straight average of the apparent column density measurements. 
In these cases the errors were determined by assuming that the continuum placement uncertainties are not 



lowered by averaging the individual measurements (we adopt the minimum), while statistical uncertainties 
are lowered by cr| /N, where Oi are the individual uncertainties and N the number of lines averaged. 

The continuum and statistical uncertainties are then added in quadrature. 

We also use the apparent column density integrations from Table 4 for deriving lower limits where 
saturation corrections seem untenable (i.e., for strongly saturated lines such as C II A1036). In these cases 
we adopt the highest (most constraining) lower limits available. We derive upper limits for undetected species 
assuming they are distributed as a Gaussian with breadths equal to those observed for other species (see 
§5.2 and §6.4). All quoted upper limits have a 3 ct confidence. We note that no molecular hydrogen is seen 
along this sight line; our limit on the H 2 column density is logA^(H 2 ) < 14.35 (3(t), summing all rotational 
states J < 3. The absence of H 2 implies the line of sight does not contain significant cold neutral gas. 

For species where saturation corrections could be important, we derive column densities by fitting a single 
component, Doppler-broadened curve of growth to the measured equivalent widths (fitting measurements 
from both detectors simultaneously) or by applying the curve of growth derived for similar species. We have 
used curve-of-growth fitting (in a manner following Howk et al. 2000) for determining column densities of 
Ar I and Fe II using all of the available measurements. We have used the resulting 6-values to provide further 
information on the degree of saturation expected for P II, S III, and Fe III. A few notes about specific species 
are in order: 

C II - The profile of the strong C II 1036.337 A transition shows the entire velocity extent of the gas 
along this sight line. We have not presented measurements of C II* at 1037.018 A for this sight line. The 
C II* absorption (which can be seen in the panel showing O VI A1037.617 in Figure 3) is contaminated by 
unidentified stellar features. While this could yield an upper limit to the C II* column density, we have not 
used the profile for this purpose given the possibility that some saturation exists in the profile (which would 
keep us from deriving a true upper limit). 

P II - The LiFlB and LiF2A measurements of P II are inconsistent with one another at the 3(7 level. The 
origins of this discrepancy are not clear, although the LiFlB velocity profile seems inconsistent with those seen 
in other species. We proceed by adopting the value derived using LiF2A, but given the discrepancies between 
the two channels, we will make little use of the P II measurements. Adopting the LiF2A measurements, 
we have derived curve-of-growth column densities assuming 5(P II) > 6(Ar I). The inequality is adopted 
because of the belief that Ar I is more likely to reside in cooler regions of the ISM. However, even adopting 
the relatively small (compared with Fe II) Ar I &-value gives results that are consistent (within ~ Icr) with 
the apparent column density determinations. For this reason we adopt the average of the apparent column 
density and curve of growth methods. 

Fe II - For Fe II we have adopted the /-values from Howk et al. (2000) with the exception of that 
for the transition at 1144.938 A, for which we have adopted the new laboratory measurement by Wiese, 
Bonvallet, and Lawler (2002). However, for comparison we have also fit a curve of growth adopting the 
FUV Fe II oscillator strengths derived by Mallouris et al. (2001). We find the column densities and 5-values 
derived with these two sets of oscillator strengths to be consistent within the Icr uncertainties. 

Fe III - We have derived Fe III column densities assuming this ion has a velocity distribution similar to 
that of Fe II, adopting 6(Fe III) = 22.9±2.0 km s“^. That is, we apply the Fe II curve of growth (with 6-value 
uncertainties twice that of Fe II) to the Fe III equivalent width measurements. We derive column densities 
in this way that are indistinguishable from those derived using the apparent column density method. 

S III - We assume that this ion is distributed in a manner like Fe III, which implies no saturation 
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correction. This is justified by the similarity of the S III and Fe III profiles (see §6.4). 

N 11 - Though this line only yields a lower limit to the N II column density, the analysis that led 
to this limit is sufficiently different than that used for the other species that it requires description. The 
strong broad stellar He II 1084.940 A transition overlaps the positive velocity wing of the N II 1083.994 A 
interstellar transition (the center of the He II is found at ~ +120 km s“^ relative to the restframe of N II, 
including the effects of the large stellar velocity). Weak absorption from stellar N II* AA1084.566,1084.584 
is also present at lower positive velocities (~ +20 km s“^) in the restframe of the ground state interstellar 
transition. We have accounted for the overlap of these stellar features with the interstellar absorption by 
using a model atmosphere provided by P. Chayer (2001, private communication) to normalize the data near 
N II A1083.994 before processing the data as described above. The model atmosphere was calculated in a 
manner similar to that described in Sonneborn et al. (2002) using the TLUSTY code of Hubeny & Lanz 
(1995) and the SYNSPEC code of I. Hubeny (2000, private communication). The model is characterized by 
T = 35,000 and logg = 4.0 (Dixon et al. 1994). The model assumes a ratio He/H = 0.1 and a metallicity 
[M/H] = -2. 


5. Highly-Ionized Gas Toward von Zeipel 1128 

The hot phase of the ISM in galaxies is a direct result of the feedback of energy and matter into the 
ISM from early-type stars. The hot (T > 10® K) ionized gas in the Galactic ISM is typically traced through 
its X-ray emission, particularly its soft X-ray emission. Much of this hot material will eventually cool 
through various processes, either through simple radiative cooling (e.g., as envisioned in galactic fountain 
scenarios; Edgar & Chevalier 1986) or through more complex mechanisms where the hot gas interacts with 
cooler material (e.g., through conduction or mixing; see review by Spitzer 1996). Material cooling from 
X-ray emitting temperatures passes very quickly through the temperature regime of a few x 10® K due to the 
presence of strong emission lines of highly-ionized metals in this temperature range, especially the lithium¬ 
like ions of carbon, nitrogen, and oxygen. For this reason, gas in this regime is often referred to as “transition 
temperature” gas or, by analogy with recent work on the intergalactic medium, warm-hot ionized gas. 

It is this transition temperature gas that is traced by O VI absorption toward vZ 1128 and other 
background ultraviolet continuum sources. This section discusses the properties of the highly-ionized gas in 
the Galactic thick disk toward vZ 1128. While O VI absorption line measurements provide good kinematic 
information on the warm-hot ionized medium, it should be emphasized that this material does not represent 
the majority of the hot ionized gas in the thick disk. Instead, O VI is a tracer of the cooling of the hot 
ionized medium or of its interaction with cooler material. 


5.1. Column Densities of Highly-Ionized Gas 

The sight line to vZ 1128 probes the entire interstellar thick disk of the Milky Way and parts of the 
more extended corona. The distribution of O VI in the former region is best described with an exponential 
scale height of hz ~ 2.3 kpc having an excess in the northern polar cap region of ~ 0.25 dex (Savage et al. 
2003; Zsargo et al. 2003). The midplane density of the distribution is no(0 VI) « 1.7 x 10“® cm“® (Jenkins, 
Bowen, & Sembach 2001). These quantities suggest an average value log N (O VI) sin |fo| of 14.08+0.25 = 14.33 
integrated through the northern extent of the Galaxy. It should be noted that the O VI layer is very patchy 
(Howk et al. 2002b; Savage et al. 2003), with a standard deviation about the mean N{0 VI) sin |5| of ~ 0.19 
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dex (~ 40%) for sight lines probing the entire Galactic halo. 

The O VI column density observed toward vZ 1128, logiV(0 VI) sin |&| « loglV(0 VI) = 14.49 ± 0.03, is 
larger than the northern Galactic average value determined by Savage et al. (2003) in their survey of O VI 
toward extragalactic continuum sources. However, given the large observed dispersion in the distribution of 
O VI column densities, that toward vZ 1128 does not seem unusually large. We note that the great distance 
of vZ 1128 makes the column densities derived here more comparable to the extragalactic sample of targets 
than, for example, to the halo star sample studied by Zsargo et al. (2003). 

In principle the column density ratios of various highly-ionized species can be used to constrain the 
temperature of the gas, if collisional ionization dominates (as we expect for O VI). For this reason we have 
searched for S VI along this sight line, which has creation and destruction ionization potentials of 72.7 and 
88.0 eV. While the O VI column density is relatively large, there is no evidence for absorption by interstellar 

5 VI in the FUSE spectrum of vZ 1128. Assuming S VI to have the same velocity dispersion as O VI, we 
derive a 3cr limit of log A^(S Yl)/N{0 VI) < — I.l. In the collisional ionization equilibrium calculations of 
solar-metallicity gas by Sutherland & Dopita (1993), this ratio implies T > 2 x 10® K. We caution that 
equilibrium conditions may not hold. For comparison, the non-equilibrium calculations of Shapiro & Moore 
(1976), in which collisionally-ionized gas at T ~ 10® K is allowed to cool isochorically to lO"* K, are consistent 
with this ratio for 10'^ < T < 6 x 10® K. Evidently the ratio N{S Yl)/N{0 VI) does not by itself place strong 
constraints on the temperature of highly-ionized gas given the discrepancies between the equilibrium and 
non-equilibrium calculations. 

Similarly, our data limit the ratio log Af(S IV)/Af(0 VI) < —0.8 (3cr). This ratio is also consistent 
with the Sutherland & Dopita (1993) equilibrium models for T > 2 x 10® K. The Shapiro & Moore (1976) 
calculations predict values in agreement with our limit for T > 1.5 x 10® K; however, their calculations for 
1.2x 10^ < T < 7.2 X lO"* K are also consistent with our limits on A^(S IV)/A^(0 VI). In this case the 
temperature limits derived from the equilibrium and non-equilibrium models overlap better than for the 
N{S Yl)/N{0 VI) ratio. However, one must still be concerned that these models are not the best descriptor 
of the situation at hand. (It is also important to note that some S IV could arise in photoionized material, 
similar to the case for the well-studied ion Si IV.) 

We have also placed limits on the column density of P V (with creation and destruction ionization 
potentials 51.4 and 65.0 eV) toward vZ 1128, deriving log A^(P Y)/N{0 VI) < —1.8 (3ct). Neither Sutherland 

6 Dopita nor Shapiro & Moore report the ionization balance of phosphorus. However, with log(P/O )0 = 
—3.29 (Grevesse & Sauval 1998), this limit is not likely to be useful. 

We can estimate the column density of protons associated with the hot ISM in this direction using the 
observed O VI column density. The proton column density associated with a tracer of ionized gas, X*, is: 

X(HII)«X(X*)(X/H)-®x(X*)”\ (3) 

where X(X*) is the column density of the ion X*, (X/H) is the abundance of X/H (usually assumed to be 
equivalent to the solar system abundance), and a;(X®) = n(X*)/n(X) is the ionization fraction of the ion X®. 
In the case at hand, we have X(H II) « X(0 VI) (0/H)~^ x{0 VI)“®. We assume the abundance of oxygen, 
0/H, to be the solar value of log(O/H )0 = —3.26 (Holweger 2001) and that the fraction of oxygen in the 
form of O VI, x(0 VI), is < 0.2, which is consistent with both the equilibrium calculations of Sutherland & 
Dopita (1993) and the non-equilibrium calculations of Shapiro & Moore (1976). With these assumptions we 
derive logX(H II) > 18.4 for the transition-temperature (~ 10® — 10® K) gas toward vZ 1128. 
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5.2. Kinematics of the Highly-Ionized Gas 

Figure 4 shows the Na{v) distribution of O VI toward vZ 1128 and the NRAO H I 21-cm spectrum. 
The middle panel shows the O VI Na{v) profile with the best-fit Gaussian distribution overlayed. The 
properties of this Gaussian are summarized in Table 6, along with the properties of Gaussians fit to the 
Na{v) profiles of several lower-ionization species. The best-fit Gaussian dispersion, cr (= bf\/2), after removal 
of the instrumental breadth, is 31.8 ± 1.1 km s“^ (equivalent to a Doppler parameter b = 45.0 km s“^). This 
limits the temperature of the O VI along the sight line to T < 2 x 10® K. This is significantly larger than 
the temperature at which the abundance of O VI peaks in collisional ionization equilibrium, T ~ 3 x 10® K, 
and the large breadth indicates significant non-thermal broadening of the profile. 

The kinematic properties of the O VI profile toward vZ 1128 are significantly different than expected for 
a highly-ionized thick disk corotating with the underlying thin disk. Because of the high latitude of this star, 
Galactic rotation is expected to have very little impact on the velocity profiles of any interstellar species. 
The bottom panel of Figure 4 shows the O VI Na{v) profile with three models of simple rotation of a thick 
O Vl-bearing disk with scale height of ~ 2.3 kpc (see Savage et al. 2003). The three models differ in the 
assumed dispersion of the O Vl-bearing gas (with models shown for Gaussian dispersion of ct = 10, 20, and 
30 km s“^). The midplane density in each model is scaled to match the red wing of the O VI profile. 

The observed velocity distribution of O VI is clearly inconsistent with these simple models of Galactic 
rotation. The observed profile shows much more material at negative velocities than is expected for a 
corotating thick disk of O Vl-bearing material. Indeed, the velocity centroid of O VI suggests the bulk of 
the highly-ionized material along this sight line is falling onto the Galactic disk from the halo. The central 
velocity of the O VI profile, (riLSu) ~ —26 km s“^, is consistent with the center of the broad negative velocity 
absorption seen in the NRAO H I profile (at wlsr ~ —25 km s“®). It is also interesting to note the relatively 
good agreement between the central velocities of WNM/WIM tracers and O VI (see Table 6). Thus, it 
would seem that the infalling gas along this sight line contains a range of ionization states. We discuss the 
kinematic relationship between the WNM/WIM tracers in §6.4. 

The kinematics of O VI along the vZ 1128 sight line are perhaps best discussed in the context of similar 
high-latitude observations. Wakker et al. (2003) and Savage et al. (2003) present Galactic O VI absorption 
lines toward ~ 100 extragalactic objects. They find that while there is a large spread in the observed central 
velocities of the O VI observed along high-latitude sight lines, the average of the distribution is ulsr = 0± 21 
km s“^ (standard deviation). Thus, while we observe the O VI to be infalling toward vZ 1128, this should 
be viewed in the context of the broad distribution of outflowing and infalling O VI gas across the entire 
high-latitude sky. The central velocity of the vZ 1128 O VI distribution cannot be used, for example, to 
claim that the entire highly-ionized Galactic thick disk is collapsing onto the plane. 

The observed Gaussian dispersion can also be discussed in the broader context of the Wakker et al. and 
Savage et al. measurements. These authors use a slightly different approach to determining the velocity 
breadths of the O VI absorption profiles, defining a quantity we refer to as b': 

^ L f(v - ll)^iVa(r^)dr^ ] 

~ [ / Na{v)dv 

These authors find (&') = 60± 15 km s“^ for their sample of extragalactic sight lines (with no correction for 
instrumental effects). Our measured dispersion corresponds to 6 « 46.5 ± 1.6 km s“® (before correcting for 
the instrumental spread function). If derived in the manner used by these authors, we find b' for the O VI 
toward vZ 1128 is 57 ± 3 km s“^. Thus, the breadth of the O VI profile toward vZ 1128 is well within one 
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standard deviation of the high-latitude mean derived by Wakker et al. and Savage et al. 


6. Warm Ionized and Neutral Gas Toward von Zeipel 1128 

One of the goals of this program is to understand the relationship between the observational tracers 
of warm neutral and ionized gas in the Galactic thick disk. This section discusses the relative elemental 
abundances in the WNM and WIM as well as the kinematics of the tracers of these phases toward vZ 1128. 

The WNM and WIM in the Galactic thick disk are probed by low and intermediate ions, respectively. 
Low-ionization species include those thought to be the dominant ionization state of each element in the 
WNM; low ions observed toward vZ 1128 include G II, N I, O I, Si II, P II, Ar I, and Fe II. We describe the 
next-higher ionization stages as intermediate ions, including the species C III, N II, S III, and Fe III observed 
toward vZ 1128. These ions are likely to arise in the WIM (and perhaps to some degree in the hot ISM; see 
Howk & Savage 1999b) along with some of the low ions. 

The low ions may arise in both the WNM and WIM (Haffner, Reynolds, & Tufte 1999; Sembach et al. 
2000 ), which is a source of some concern, since these are the ions used to determine gas-phase abundances 
in the WNM (see discussion in Sembach et al. 2000). The degree of uncertainty associated with this 
“contamination” (or confusion) of WNM species by gas associated with the WIM varies from element to 
element. We discuss such contamination for the specific case of iron below. 


6.1. Gas-phase Abundances in the Warm Neutral Material 

The gas-phase abundance of an element X relative to an element Y is calculated as 

[X/Y]=\og[NiX^)/NiY^)]-log{X/Y)Q-log[xiX^)/x(Y^)]. (4) 

In this expression A(A*) and N{Y^) are the column densities of the ions A* and YY These column densities 
are normalized to the relative solar system abundances of the two elements, {X/Y)q (Grevesse & Sauval 
1998; see Table 7). The quantities x{X'^) and x{Y^) are the ionization fractions of each ion in the gas being 
studied. The last term is the ionization correction. Typically the dominant WNM ionization state of the 
two elements is used, and the ionization correction is assumed to be insignificant. Exceptions to this general 
rule include studies of Ar I, which is easily photoionized in WNM conditions (Sofia & Jenkins 1998). 

Table 7 and Figure 5 summarize the gas-phase abundances in the WNM derived in the usual way. 
That is, we have compared the column densities of low-ionization species with that of H I in this direction. 
Studies of the WNM typically use H I or a non-depleted metal ion as a reference for studying interstellar 
grain composition. The relative deficit of metals in the ISM compared with the solar system is assumed to 
be caused by the incorporation of these metals into dust grains (the solid phase). Applying the standard 
technique yields only an approximate determination of [A/H] for the WNM in this direction because of the 
assumption that log [a;(A*)/a;(H I)] = 0 in Eqn. (4). 

Figure 5 shows the WNM gas-phase abundances for the vZ 1128 sight line along with the range of values 
seen along other thick disk (or “halo”) sight lines. The gas-phase abundances toward vZ 1128 are not out 
of the ordinary for high-latitude sight lines. We point out that these measurements represent averages over 
the entire sight line. There is evidence that many of the metals are distributed differently in velocity space 
than the H I, especially the species heavily incorporated into dust grains (see below). 
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6.2. Ionization Effects 

The measurements given in Table 7 provide a baseline for determining the extent to which ionization 
effects influence such abundance determinations and for comparing the properties of the WNM and WIM 
along the vZ 1128 sight line. We first discuss the average ionization state of the gas along this sight line 
because it bears on the gas-phase abundances derived above. 

Using the Eqn. (3), it is possible to estimate the column density of protons associated with the interme¬ 
diate ions toward vZ 1128. Such an estimate, combined with measures of A'^(H I), provides a rough measure 
of the total hydrogen column density along the sight line. For determining the proton column density of the 
WIM, we make use of the S III column density reported in Table 5. 

The most basic assumption, i.e., that all of the WIM sulfur is in the form of S III [a;(S III) = 1] and sulfur 
is undepleted (Howk, Savage, & Fabian 1999), results in a lower limit to the proton column associated with 
the WIM: log A^(H II) > 19.1. Sembach et al. (2000) have presented models of the ionization of the WIM 
assuming ionization by early-type stars. If we adopt the sulfur ionization fractions from their composite 
standard model [a:(S III) ~ 0.18] for our calculation, we find logfV(HII) = 19.9. For comparison, similar 
calculations performed using the lower and upper limits provided by N II and P III yield 18.6 < logfV(H II) < 
20.8, consistent with the valued derived using S III. 

The proton columns estimated using S III above imply a total hydrogen column density of log fV(H U-H II) 
20.0 to 20.2 and (a;(H II)) = fV(H II)/7V(H U-H II) > 0.12 to ~ 0.46, with the latter being the preferred es¬ 
timate. These values are roughly consistent with the estimates of a:(H II) derived from the Ha intensity in 
this general direction (§2). 

With this estimate of the total hydrogen column density of the WNM and WIM (since the hot ISM 
likely contributes little, as discussed in §5), we can further estimate the total gas-phase iron abundance. 
The total iron column density along this sight line is log N(Fe)totai ~ logN(Fe lU-Fe III) = 14.95 ± 0.04, 
assuming higher ionization states represent a small fraction of the total. Thus, the gas-phase abundance of 
iron along the sight line to vZ 1128 (averaged over all phases and regions) is 

-0.76 < [Fe/H],„,^, < -0.56. 

This can be compared to the value [Fe/H]wArM = —0.64±0.06 derived using Fe II and H I with no ionization 
correction. The ionization corrections for the derived WNM gas-phase abundance appear to be < ±0.1 dex 
(averaged over the entire sight line). 


6.3. Gas-phase Abundances in the Warm Ionized Material 

It is also possible to estimate the gas-phase abundances of refractory elements in the WIM by comparing 
the columns of intermediate-ionization stages of depleted and non-depleted species and applying an ionization 
correction. Howk & Savage (1999b) have presented determinations of [Al/S] in the WIM using this approach, 
arguing on the basis of the derived values that dust is present in the ionized medium of the Galaxy. This 
conclusion has been confirmed by the detection of far-infrared emission from dust associated with the WIM 
(Lagache et al. 2000). 

Our observations of Fe III and S III provide a similar opportunity to determine the relative gas-phase 
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abundances of iron and sulfur in the WIM. We determine [Fe/SJiy/M in the manner outlined in Eqn. (4): 

[Fe/S']wiM = log [N(Fe III)/iV(S III)] - log (Fe/S)© - log [a;(Fe III)/x(S III)]. (5) 

In this case the ionization correction term, log]x(Fe III)/x(S III)], is significant. 

The composite standard WIM model (Sembach et al. 2000) has log[a;(Fe III)/a;(S III)] = +0.45, which 
yields 

-0-74 ±0.3 

for the sight line to vZ 1128. The error bar is an estimate of the uncertainties in the ionization corrections due 
to the sensitivity of a;(Fe III)/a;(S III) to the adopted effective temperature of the ionizing source following 
Howk & Savage (1999b). This determination is in excellent agreement with the WNM determination and 
the values of [Fe/Hj^^^^j and derived above. The implications of this similarity are discussed 

more thoroughly in §8. 


6.4. Kinematics of the Warm Ionized and Warm Neutral Gas 

We now turn our attention to the kinematic relationship between the low- and intermediate-ions toward 
vZ 1128. Figure 6 shows the apparent column density profiles of a range of ions observed toward vZ 1128, 
including species tracing the WNM and the WIM. Also shown is the O VI profile, which traces the warm-hot 
ISM. Measurements of the velocity centroids and breadths of the Na{v) profiles of several species are given 
in Table 6. The peak apparent optical depths of each transition are also given as a guide for judging the 
relative importance of saturation effects on the derived kinematic properties. 

It is clear that the resolution of FUSE is insufficient to clearly resolve the detailed component structure of 
this sight line. Examining the Na{v) profile of the strong Si II 1020 A line (which is susceptible to unresolved 
saturation at the highest optical depths), there is clear evidence for velocity substructure that cannot be 
clearly separated at FUSE resolution. For this reason, we concentrate our discussion of the kinematics of 
the WNM/WIM gas on the central velocity and breadth of the absorption. The majority of the species 
presented in Figure 6 and Table 6 have similar central velocities ranging from ulsr ~ ~31 to —21 km s“^, 
including O VI. The one exception to this agreement in the velocities of the observed species is Ar I, which 
is centered near vlsr ~ —10 km s“^. Because of its large ionization cross section (Sofia & Jenkins 1998), 
Ar I is relatively easy to ionize in diffuse interstellar gas, even though its ionization potential is larger than 
that of H I. This is the most likely cause for the apparent deficit of Ar I seen in Figure 5, particularly since 
argon is unlikely to bind to grains effectively (Sofia & Jenkins 1998). All of the species displayed in Figure 
6 show absorption at the velocities at which Ar I has a maximum optical depth. 

To the extent that the FUSE resolution allows us to study such coincidences, the breadths of the profiles 
of species tracing the WNM and WIM are quite similar. Figure 7 compares the apparent column density 
profiles for Fe II, Fe III, and S III. With the exception of a possible weak excess in Fe III at wlsr ~ —50 
km s“^, the profiles of these species track one another remarkably well. The similarity of the observed 
profiles suggests the WNM and WIM share a similar component structure. Large component-to-component 
variations in the ratio of the WNM and WIM tracers can be ruled out for the dominant component groups 
along this sight line (strong variations in the ratio among low column density components cannot be ruled 
out). 

An interpretation of the kinematics of the WNM/WIM tracers discussed in this section is that the species 
centered near ulsr ~ — 26 km s“^ principally trace gas in the Galactic thick disk and are associated with the 
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broad H I seen at similar velocities. The material closer to the LSR zero point, including the narrower peak 
of the H I distribution and the majority of the Ar I, is likely tracing gas associated with lower-ionization, 
higher-depletion material in the Galactic disk in this direction. The principal WNM tracers accessible by 
FUSE probe elements that are readily incorporated into dust grains. These ions are therefore biased toward 
gas in the Galactic thick disk, which is more likely to be ionized than the thin disk and to have higher 
gas-phase abundances due to the processing of dust grains (e.g., Sembach & Savage 1996). 

It is signihcant that there seems to be little evidence for gas that can be associated with only the neutral 
or ionized gas of the Galactic thick disk. 


7. The Lack of High-Velocity Gas 
7.1. High Velocity Clouds 

There is no evidence for discrete high-velocity cloud material along the vZ 1128 sight line in O VI or 
lower-ionization species. Sembach et al. (2003) have discussed the presence of high-velocity O VI toward 
extragalactic objects. They find high-velocity O VI toward ~ 60% of the sight lines studied (although some 
were initially chosen to probe the known H I high-velocity clouds). In a few cases, the lack of high-velocity 
material toward vZ 1128 is confused by the presence of high negative-velocity stellar photospheric absorption 
features (see, e.g., the C II 1036.337 A transition in Figure 3). 

For O VI, our 3cr upper limit on the column density of high-velocity (|wlsrI ^ 125 km s“^) material is 
log V(0 VI) < 13.3, assuming cr ~ 32 km s“^ {h ~ 45 km s“^). We note that ~ 20% of the sight lines studied 
by Sembach et al. (2003) show the presence of very broad, high-positive-velocity wings in O VI extending 
from |wlsr| ~ 100 to ~ 250 km s“^. Most of these wings are found in the northern Galactic sky. We see no 
evidence for low optical depth, high-velocity wings such as these toward vZ 1128. 

There are no H I high-velocity cloud complexes toward vZ 1128 identified through 21-cm H I emission 
line observations. Thus, even though a large fraction of the extragalactic targets show high-velocity O VI, 
it is unclear if we expect to see such material toward vZ 1128. To demonstrate the lack of high-velocity 
material toward vZ 1128, we show normalized spectra of two higher-order Lyman series H I lines toward 
vZ 1128, together with the proHles of G II A1036.337, G III A977.020, and O VI A1031.926 in Figure 8. The 
H I proHles have been normalized by the same model stellar atmosphere used for the N II profile (see §4). 
Even though this model is a reasonable descriptor of the broad stellar H I absorption seen in the FUSE data, 
we caution that the application of an incorrect model profile could leave unrealistic broad absorption in the 
spectrum, particularly at negative velocities where the absorption from the star is strongest. The vertical 
lines in Figure 8 denote the approximate range of C II absorption at its half intensity points. 

There is little visible H I beyond the extent of the C II proHle. The H I proHles seem to be confined 
to —125 < ulsr ^ +65 km s“^ (with the caveat that the negative velocity extent is dependent upon the 
adopted stellar profile). Outside of this velocity range, we find an upper limit to the H I column density of 
any high-velocity material of log A^(H I) < 14.7 (3(j). Neither do we see evidence for high-velocity absorption 
in low-ionization metal species. Assuming such absorption would have breadths similar to the values listed 
in Table 6, we derive limits of logA^(Fe II) < 13.4, log A^(Fe III) < 13.5, and logA^(N II) < 13.4 (Sa limits) 
for material with [ulsrI ^ 125 km s“^. The Fe III limits only apply to high positive velocities since the 
profile is blended with interstellar Fe II and stellar Fe III at negative velocities. 

The intermediate-velocity cloud complex K (Wakker 2002) is found at lower Galactic latitudes in H I 
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maps at t^LSR ~ ~70 km s“^. The map of this complex in Wakker (2002) shows the nearest contour enclosing 
A^(H I) > 5 X 10^® is ~ 10° from vZ 1128. It is difficult to say if complex K is detected toward vZ 1128. The 
very strongly-saturated C II 1036.337 A transition shows some evidence for weak absorption at velocities 
consistent with an origin in complex K. Small wings at the appropriate velocity may also be seen in the 
stronger O I and Si II transitions (see Figure 3). However, the absorption at these velocities is not sufficiently 
well delineated to allow us to study it in detail. 


7.2. Globular Cluster Gas 

We note that no absorption associated with the circumstellar environment of vZ 1128 or with the 
globular cluster M 3 is seen in our data, assuming any gas in the vicinity of M 3 would be found near its 
systemic velocity (ulsr = —137.0 ± 0.3; Soderberg et al. 1999). For O VI, our 3a upper limit on the column 
density of material associated with the globular cluster is log7V(0 VI) < 13.3, assuming cr ~ 32 km s“^ 
{b ~ 45 km s“^). Neither do we see evidence for absorption in lower-ionization species. Assuming such 
absorption would have breadths similar to the values listed in Table 6, we derive limits of log A^(Fe II) < 13.4 
and logA^(N II) < 13.4 (3 ct). The Lyman series of hydrogen cannot be used to search for globular cluster or 
circumstellar material because of the strong stellar H I absorption. 

These determinations can be used to estimate upper limits on the mass of interstellar material in the 
globular cluster. The calculation follows that for determining the hydrogen column density associated with 
O VI in the hot ISM (see §5) with the additional step of multiplying that column by a projected area for 
the cluster. We assume a metallicity [M/H] « —1.57 (Harris 1996) and that the upper limits of the above 
ions are applicable to the whole cluster. The distance of vZ 1128 from the cluster core (~ 3.'9 or ~ 11 
pc) and the limits on O VI imply an intracluster gas mass (M/cm) of M/cm ^ 4.4a;(0 VI)“^ Mq. If we 
limit this calculation to temperatures 5.3 < logT < 5.6 [the range over which a;(0 VI) > 0.05 for both the 
equilibrium calculations of Sutherland & Dopita (1993) and the non-equilibrium calculations of Shapiro & 
Moore (1976)], we have M/cm ^ 90 M©. Similarly we can use the limits on N II assuming the gas is fully 
ionized (i.e., hydrogen is fully ionized and the ionization of nitrogen is tied to that of hydrogen via its charge 
exchange reaction) to limit the mass of the intracluster medium to be M/cm ^ 35 M©. 

Our calculations assume that the limits to the gas column density toward vZ 1128 are applicable to the 
entire projected area within the 11 pc distance to the cluster core.® If the column of gas is higher toward the 
center of the cluster or if vZ 1128 lies on the front side of the cluster, then our limits are not appropriate. 
Our limits on the mass of ICM material in M 3 are similar to other recent limits on the amount of ionized 
gas in globular clusters (e.g., Knapp et al. 1996) and the one detection (Freire et al. 2001). The density of 
ionized gas detected in 47 Tuc by Freire et al. (2002) implies a mass of ~ 10 Mq if the gas extends as far as 
11 pc from the cluster center, the projected distance of vZ 1128 from the center of M 3. 


7.3. High-Velocity Dispersion Gas 

Kalberla et al. (1998) have used the LDS H I observations to argue for the presence of a high velocity- 
dispersion component of the Galactic ISM. They find evidence for diffuse H I with cr « 60 km s“^ and 


®We note that M 3 has a half-mass radius of ~ l.'l or ~ 3.3 pc and a tidal radius of ~ 0.'5 or ~ 1.6 pc (Harris 1996). 
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logA^(HI) « 19.15. If exponentially distributed, Kalberla et al. argue such a distribution would have a 
midplane density of ~ 1.2 x 10“^ cm“^ and a scale height of ~ 4.4 kpc. 

Such a component of the ISM should be readily visible as strongly-saturated Lyman-series absorption 
toward distant stars and AGN. Savage et al. (2000) have argued on the basis of the measured equivalent 
widths of the strong Mg II doublet near 2800 A that such a phase of the ISM cannot exist unless the 
abundance of Mg in this material is far below solar (~ 0.05 solar). 

Using the Lyman series of H I as a probe of this gas alleviates the issue of abundances. If such material 
were present toward vZ 1128, the H I 930.748 A line seen in Figure 8 would have < 1% residual intensity 
over the range —200 < 'Clsr ^ +200 km s“^. The H I Lyman series towards vZ 1128 (conservatively) limits 
the column density of such material in this direction to be log Af(H I) < 16.7 cm“^. 

While we can say with some confidence that a phase of the ISM which is smoothly distributed and 
has a dispersion and column density like those discussed in Kalberla et al. (1998) does not exist toward 
vZ 1128, their detection as a whole cannot be ruled out by the present data alone. To identify this phase 
of the ISM, Kalberla et al. averaged over ~ 10° x 10° sections of the sky to identify this gas. If this high 
velocity-disperion H I distribution is contained in clouds which cover a relatively small fraction of the sky and 
as an ensemble have a cloud-to-cloud dispersion of order 60 km s“^, we would not necessarily expect to find 
material associated with this distribution along a single sight line. Lockman (2002) has recently idenfitied 
a population of discrete thick disk clouds in the inner Galaxy. Whether the ensemble of such clouds at the 
solar circle (if such a population exists) can match the distribution inferred by Kalberla et al. (1998) is not 
known. 

An investigation of the validity of the Kalberla et al. (1998) result is best pursued using a large 
database of FUSE observations of extragalactic sources (where there is no broad Lyman-series absorption in 
the background source itself). Should no such sight lines through the entire Galactic halo show evidence for 
such gas, the reality of such an H I layer will need to be reevaluated. 


8. Discussion 

8.1. The Phases of the Galactic Thick Disk 

Because we are able to measure the column densities of species associated with all phases of the ISM 
toward vZ 1128, we can estimate the contribution each phase makes to the total hydrogen (neutral plus 
ionized) along this sight line. In this way we can construct a census of the phases in this direction. Table 
8 gives such a summary of the phases toward vZ 1128, which contains the hydrogen column for each phase 
(as estimated above) and the fraction of the total each represents along this sight line. We also summarize 
the best estimates of the scale heights of each phase from the literature. 

The fraction of the total hydrogen along the vZ 1128 sight line estimated to be associated with warm 
ionized gas is ~ 10% — 45%. In this case, the higher end of this range is more robust since it is based upon a 
model of the ionization fraction of sulfur in the WIM. The lower value assumes that all of the sulfur is in S III, 
which yields a firm lower limit to the ionized hydrogen column density associated with the WIM. This range 
of fractional hydrogen column density is similar to the observed values derived by comparing the electron 
column densities (derived from pulsar dispersion measures) with H I column densities toward distant globular 
clusters. Reynolds (1991b) summarizes five such measurements which give fractions of ionized hydrogen in 
the range 19% — 40%. 



- 18 - 


The highly-ionized, transition-temperature gas traced by O VI represents a relatively small fraction of 
the total hydrogen column density along this path through the Galactic thick disk and corona. It should 
be noted that the O VI absorption is likely not tracing material at X-ray emitting temperatures (> 10® K); 
therefore, the total amount of hot gas could be higher. Furthermore, because we have adopted the highest 
ionization fraction for estimating the hydrogen column from O VI, the stated values should be treated as 
(uncertain) lower limits. 


8.2. The Relationship Between Neutral and Ionized Gas in the Galactic Thick Disk 

Our determinations of the depleted gas-phase abundances of iron in the WNM and WIM are very similar, 
suggesting the content of iron-bearing dust is similar in these two phases. This tends to support the notion 
that the gas and dust that make up the WNM and WIM of the Galactic thick disk have experienced similar 
dust destruction and formation histories. If the processes that feed WNM and WIM material into the thick 
disk were significantly different (i.e., if one were more violent than the other), the dust content, and hence 
the gas-phase abundances, of these phases might be significantly different. 

We point out, however, that our determination of has a relatively large uncertainty due to 

the photoionization modeling. Using refined models and species that are less susceptible to changes in the 
ionizing spectrum (e.g., A1 III; see Howk & Savage 1999b) for determining the WIM abundances will provide 
more precise results. It is also worth pointing out that the WNM abundance determinations may be affected 
at the ~ 0.1 dex level by our neglect of the ionization correction term in Eqn. (4). 

The kinematic similarities between the tracers of the WNM (e.g., Fe II) and the WIM (e.g., Fe III, S III) 
along this sight line also suggest a relatively close relationship between these phases of the Galactic thick 
disk. The distribution of clouds in velocity space is likely very similar between the two phases. We see no 
significant difference in the large-scale velocity dispersion of the warm ionized gas and warm neutral gas 
toward vZ 1128. This indicates that these two phases have similar vertical extents along this high-latitude 
sight line, unless processes are at work that affect the vertical distribution of two phases differently while 
suppressing large-scale velocity profile differences. 

The kinematic similarities suggest a relationship between these phases toward vZ 1128 that is similar 
to that seen toward HD 93521 (Spitzer & Fitzpatrick 1993) and p Leo (Howk & Savage 1999b). Along these 
sight lines, tracers of electron-rich gas (e.g., G II*, A1 III, S III) show a velocity component structure that 
is virtually indistinguishable from tracers of the WNM (e.g., Zn II, Si II, Fe II). The apparent similarity of 
the Fe II and Fe III profiles as observed by FUSE does not rule out moderate, unresolved component-to- 
component variations in the ratio of these species. 

The apparent constancy of the Fe III to Fe II ratio with velocity also implies that there are no strong 
intermediate- or high-velocity features with unusual Fe III content (i.e., none that are significantly more 
highly-ionized than the bulk of the thick disk material along this sight line). The resolution of FUSE does 
not allow us to place strong limits on such components at low velocities, and the overlap of the strong Fe II 
II2I.975 A transition at negative velocities relative to Fe III also makes the limits more significant at positive 
velocities. 

There is one caveat to these kinematic comparisons of the WNM and WIM tracers. If the WIM does 
indeed represent as much as ~ 45% of the total hydrogen in this direction, much of the column density of 
singly-ionized species we have associated with the WNM may arise in the WIM. For example, the standard 
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composite WIM model of Sembach et al. (2000) predicts a;(Fe II) « 0.5. If the fraction of hydrogen toward 
vZ 1128 is ~ 50%, then ~ 25% of the Fe II is expected to arise in the WIM. In the case of Si II, a;(Si II) ~ 0.9, 
implying that 40% of the column density of this ion could arise in the WIM. It is not difficult to understand 
why the tracers we have naively associated with the WNM might have very similar kinematics to the 
intermediate ions if they are indeed contaminated by this amount of WIM gas. 

The observed O VI profile, which traces transition-temperature gas in the warm-hot ISM, is centered at 
velocities consistent with the centroids of the WNM/WIM tracers, but has a much larger velocity dispersion 
compared with the tracers of the low and intermediate ions. The large breadth of the O VI profile may 
be the result of large turbulent (or other non-thermal microphysical) motions or large-scale cloud-to-cloud 
motions. The breadth of the O VI profile toward vZ 1128 is consistent with the N{0 VI)-line width relation 
identified by Heckman et al. (2002). 

The agreement between the central velocities of O VI and the lower ionization species is intriguing, 
though perhaps coincidental. The correspondence of central velocities of all of the tracers of the thick disk 
ISM in this direction is naturally explained if the O VI arises in interfaces between the warm gas (neutral and 
ionized) and a much hotter medium (X-ray emitting material at T > 10® K). Because of the short cooling 
time of gas in the range of temperatures at which O VI is abundant, many classes of theoretical models 
invoke such an arrangement (see summary in Spitzer 1996). [Indeed, such an arrangement calls to mind the 
clouds envisioned by McKee & Ostriker (1977).] 

Examining Figure 8 one finds that the total velocity extent of the strong C II 1036.337 A and O VI 
1031.926 A profiles is similar. That is, most of the absorption seen in O VI is seen at velocities where 
there is also some lower-ionization material (as traced by C II). The much larger apparent dispersion of the 
O VI profile compared with lower-ionization lines of similar peak optical depths is likely due to components 
further from the velocity centroid which exhibit higher O VI to low ion column density ratios. If one assumes 
a relatively constant quantity of O VI per interface (e.g., Shelton & Cox 1994), then the observed O VI 
profile traces the distribution of interfaces as a function of velocity. The individual interfaces much overlap 
in velocity so that they are not visible as individual components (perhaps being within Av ~ 2a of one 
another, or ~ 25 km s“^ for T ~ 3 x 10®). 

Howk et al. (2002b; see also Howk et al. 2002a) have noted that the distribution of O VI and Fe II 
absorption toward stars in the Large Magellanic Cloud [{I, h) « (280°, —33°)] also requires that the O VI be 
closely associated with low-ionization material. In these directions, the centroids of the Milky Way O VI 
profiles are shifted relative to the Fe II profiles, although the total extent of the O VI and Fe II absorption 
seem to be well correlated. That is, the O VI and Fe II absorption seemed to trace many of the same 
structures, though the changing O VI/Fe II ratio causes the velocity centroids to be different. 

If the O VI and low-ionization species trace different portions of the same structures, the larger dispersion 
of O VI cannot be simply explained by a larger temperature in the O Vl-bearing region. Significant variations 
must exist in the low-ion to O VI column density ratio, with the outlying O VI components tracing higher- 
ionization gas (i.e., material with a higher O VI to low ion ratio). This is required not only so that the 
velocity dispersion be lower in the WNM/WIM tracers, but also to be consistent with the larger O VI 
scale height compared with the WNM or WIM. Given the much larger scale height of the O Vl-bearing gas 
compared with the WIM [hz{0 VI) « 2.3h2(WIM); see Table 8], the ratio of O VI to WNM/WIM tracers 
cannot be constant with height; a larger ratio is required for material further from the plane. 

Furthermore, if one solely invokes differences in the temperature of the O Vl-bearing gas, the required 
temperatures are too high to maintain large amounts of O VI. The largest instrumentally-corrected dispersion 
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in the profiles of lower-ionization species with peak optical depths similar to that of O VI (in order to avoid 
complications from strong saturation) is that of S III: (Tq = 14.7± 1.7. If the observed O VI is associated with 
the same components making up the S III profile, an additional broadening of cr ~ 28 km s“^ is required to 
explain the breadth of the O VI profile. This would require T ~ 1.5 x 10® K in each component. At such 
high temperatures O VI is not expected to be abundant (Sutherland & Dopita 1993; Shapiro & Moore 1976). 

The observed arrangement of interstellar absorption toward vZ 1128 can perhaps be best explained by 
a scenario in which O VI and the lower-ionization species arise in the same structures, with the O VI arising 
in an interface between the warm and hot phases of the ISM. The material further from the line centroid is 
found further from the plane of the Galaxy and is characterized by a higher ionization state (as one would 
expect for lower warm material per constant-column O VI interface). In this scenario, the coincidence of the 
velocity centroids of O VI and lower-ions requires that most WNM/WIM clouds of the Galactic thick disk 
have a corresponding O Vl-bearing interace, implying that each cloud is embedded within a hot medium. 
Further analysis of the relationship between O VI and the lower-ionization gas along a larger sample of sight 
lines is required to determine whether or not this agreement is simply coincidence. 


8.3. Ionization Corrections to Gas-Phase Abundance Determinations 

As discussed in §6.1, the [Fe/H] abundances derived using the traditional measures (i.e., Fe II/H I) 
provide a determination that is within ~ ±0.1 dex of the value derived by accounting for all of the pertinent 
ionization stages of iron and hydrogen [i.e., (Fe II ± Fe III)/(H I ± H II)]. This allowable range of ionization 
corrections for determinations of [Fe/H] is relatively small, implying little ambiguity in studies of gas- 
phase abundances in varying Galactic environments. For iron, which is very heavily incorporated into dust 
grains, this ionization correction is indeed small compared with the typical effects of dust. Iron is typically 
underabundant (depleted) by more than 0.5 dex for WNM material in the Galactic halo and 1.0 dex for 
WNM material in the Galactic disk (Savage & Sembach 1996). However, it is important to note that the 
uncertainties due to the neglect of ionization corrections are much larger than typical error bars for such 
measurements. 

The ionization corrections implied here for studying iron in the WNM agrees well with the limits to the 
ionization corrections derived by Sembach & Savage (1996). We also note reasonable agreement with the 
ionization corrections estimated for the sight line to HD 93521 by Sembach et al. (2000), depending on the 
specific model adopted (see their Table 7). However, we caution that the Sembach et al. discussion of the 
HD 93521 sight line applies to corrections for the ions of sulfur; it is not necessarily the case that ionization 
corrections derived for iron will apply to the full range of observable singly-ionized species when compared 
with H I. In particular, the ionization fraction of Fe II in the WIM is expected to be significantly lower 
than other singly-ionized species such as S II, Si II, and P II (Sembach et al. 2000). Thus, the ionization 
corrections for these species may be larger than that for iron. 

More extensive studies of the ionization corrections along high-latitude sight lines would clearly be useful 
in determining the applicability of our values to broader studies of the Galactic thick disk. 
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9. Summary 

We have presented Far Ultraviolet Spectroscopic Explorer observations of the distant post-AGB star 
von Zeipel 1128, which resides in the globular cluster Messier 3. We have derived column densities and 
kinematic information for the prominent phases of the ISM along this sight line using the absorption lines 
present in the FUSE data. The major conclusions of this work are as follows. 

1. Thick disk material toward vZ 1128 is centered near wlsr ~ ~26 km s“^, implying infall of this material 
onto the Galactic plane in this general direction. Species tracing the prominent warm neutral, warm 
ionized, and transition temperature ionized phases are seen at very similar velocities along this sight 
line. 

2. As much as ~ 45% of the hydrogen toward vZ 1128 may be associated with the WIM. The tracers 
of WNM and WIM gas toward vZ 1128 have very similar kinematic profiles, showing both the same 
central velocity and very similar breadths. 

3. After accounting for all of the ionization stages of iron and hydrogen, we find that the ionization 
corrections to the standard method of determining the gas-phase abundance of [Fe/H] are likely of 
order 0.1 dex. 

4. The gas-phase abundance of [Fe/S] in the WIM is estimated by comparing the ratio Fe III/S III and 
applying an ionization correction (using the models of Sembach et al. 2000). The gas-phase abundance 
of iron in the WIM is indistinguishable from that of the WNM toward vZ 1128. The WIM gas-phase 
abundance of iron, [Fe/S]w/M ~ —0.74 ± 0.3, suggests significant incorporation of iron into dust in 
this phase of the Galactic thick disk. 

5. The similarities of the kinematics and gas-phase abundances of the WNM and the WIM suggest the 
two phases of the Galactic thick disk are closely related along this sight line. The WNM and WIM in 
this direction likely have similar vertical extents and similar Fe-bearing dust content. 

6. The broad (cr = 32 km s“^) O VI absorption toward vZ 1128 is centered at the same velocities as the 
WNM and WIM tracers along this sight line. The low- and high-ionization gas along this sight line 
seem to be closely related given the similarity in their central velocities and total velocity extent. It 
seems likely that the O VI arises in interfaces between hot and warm gas, with the shape of the O VI 
profile tracing the distribution of interfaces in velocity. 

7. We see no high-velocity material that could be associated with the circumstellar environment of vZ 1128 
or with the globular cluster Messier 3 in which vZ 1128 resides. No high-velocity clouds are seen along 
this sight line in the low- or high-ionization states observed by FUSE. No high-velocity dispersion gas 
is seen in the H I Lyman-series absorption. 


We thank P. Ghayer for calculating a model stellar atmosphere for comparison with the EUSE obser¬ 
vations of vZ 1128. We make use of data from the Wisconsin H-Alpha Mapper, which is funded by the 
National Science Foundation. This work is based on data obtained for the Guaranteed Time Team by the 
NASA-CNES-GSA FUSE mission operated by the Johns Hopkins University. Financial support to U. S. par¬ 
ticipants has been provided by NASA contract NAS5-32985. JGH and KRS recognize support from NASA 
Long Term Space Astrophysics grant NAG5-3485 through The Johns Hopkins University. 
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Table 1. Properties of von Zeipel 1128 


Property 

Value 

Reference®' 

Cluster 

M 3 (NGC 5272) 


RA (J2000) 

13 42 16.84 


Dec. (J2000) 

-k28 26 00.8 


{l,b) 

(42?5, -f78?7) 


Distance [kpc] 

10.2 

1 

2 [kpc] 

10.0 

1 

Sp. Type 

08p 

2 

V [mag.] 

14.93 

3 

B-V [mag.] 

-0.27 

3 

Teff [K] 

35,000 ± 1,000 

4 

logff 

4.0 ±0.25 

4 

log L/Lq 

3.21 ±0.12 

4 

ULSR, cluster [km S“^] 

-137.0 ±0.3 

5 

ULSR,vZ1128 [km S“\ 

[ -140 ±8 

6 


^'References: (1) Djorgovski 1993; (2) Garrison & Al¬ 
bert 1986; (3) Johnson & Sandage 1956; (4) Dixon et 
al. 1994; (5) Soderberg et al. 1999; (6) this work. 







Table 2. Log of FUSE Observations 


Visit 

ID 

Date 

No. of 
Exp. 

Exp. Time 
[sec] 

P1014101 

6/18/2000 

7 

9119 

P1014102 

6/19/2000 

9 

7684 

P1014103 

6/22/2000 

9 

14723 


Total: 

25 

31526 


Table 3. Stellar Radial Velocity Measurements^ 


Species 

A 

Ref.'" 

ULSR [km s ''j 

P IV 

1030.515 

1 

-140.7± 2.3 

P IV 

1033.112 

1 

-135 ±7 

P IV 

1035.516 

1 

-136 ±3 

SIV 

1062.678 

2 

-142.1 ±1.1 

SIV 

1072.996 

2 

-138.3 ± 1.1 

SIV 

1073.528 

2 

-139.1 ±1.0 

P V 

1117.977 

3 

-141.6 ± 1.5 

P V 

1128.008 

3 

-141.6 ± 1.3 

C III 

1174.933 

4 

-139.2 ± 1.4 

C III 

1176.370 

4 

-141.1 ± 1.3 


'^We have chosen only cleanly-separated 
transitions in the LiFlA and LiFlB spectra 
for these measurements. 

Wavelength references: (1) Zetterberg & 
Magnusson 1977; (2) Irwin & Livingston 
1976; (3) Magnusson & Zetterberg 1974; (4) 
Moore 1970. 
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Table 4. Selected Interstellar Absorption Lines Towards vZ 1128®' 


Species 

A 

[A] 

log A/ 

ITa 

[mA] 

log Na (v) 

V-,V+ 

[km s"'] 

[Det 1] 

[Det 2] 

[Det 1] 

[Det 2] 

CII 

1036.337 

2.106 

399 ±5 

400 ±7 

> 14.55 

> 14.55 

-120, ±45 

cm 

977.020 

2.872 

457 ±9 

452 ±7 

> 13.85 

> 13.85 

-100, ±120 

N I 

964.626 

0.959 

110 ±6 

115 ±5 

15.28 ±0.03 

15.35 ±0.03 

-65, ±45 

N I 

1134.165 

1.238 

154 ± 4'’ 

172 ± 3*^ 

15.22 ±0.03 

15.33 ±0.03 

-55, ±50 

N I 

1134.415 

1.528 

207 ±5 

204 ±4 

15.06 ±0.03 

15.15 ±0.04 

-50, ±50 

N I 

1134.980 

1.693 

220 ±6 

231 ±4 

14.94 ±0.04 

15.08 ±0.18 

-95, ±45 

N IF 

1083.994 

2.100 

393 ± 11 

393 ± 13 

> 14.52 

> 14.52 

-110, ±55 

O I 

924.950 

0.155 

134 ±9 

145 ± 7 

16.25 ±0.04 

16.37 ±0.03 

-65, ±40 

O I 

929.517 

0.329 

150 ± 12 

167 ±7 

16.16 ±0.04 

16.25 ±0.03 

-95, ±50 

O I 

950.885 

0.176 

146 ±8 

157 ±6 

16.23 ±0.03 

16.34 ±0.02 

-65, ±30 

O I 

1039.230 

0.980 

249 ±4 

245 ±5 

> 15.89 

> 15.84 

-85, ±40 

O VI 

1031.926 

2.137 

260 ±7 

263 ±8 

14.49 ±0.03 

14.51 ±0.03 

-160, ±100 

O VI 

1037.617 

1.836 

150 ±6 

150 ±7 

14.48 ±0.03 

14.47 ± 0.03 

-95, ±50 

Si II 

1020.699 

1.460 

153 ±4 

158 ±6 

15.20 ±0.02 

15.22 ±0.03 

-80, ±30 

P II 

1152.818 

2.435 

38 ±6'' 

60 ±4'' 

13.20 ±0.07'' 

13.41 ± 0.04'' 

-75, ±30 

P III 

998.000 

2.047 

< 34 

< 25 

< 13.53 

< 13.41 


P V 

1117.977 

2.723 

< 25 

< 25 

< 12.68 

< 12.68 


S III 

1012.495 

1.556 

92 ±4 

91 ±6 

14.47 ±0.02 

14.47 ± 0.03 

-90, ±20 

S IV 

1062.664 

1.799 

< 28 

< 34 

< 13.67 

< 13.75 


S VI 

944.523 

2.311 

< 43 

< 45 

< 13.40 

< 13.40 


Arl 

1048.220 

2.408 

90 ±4 

97 ±4 

13.70 ±0.03 

13.71 ±0.02 

-70, ±25 

Arl 

1066.660 

1.851 

51 ±3 

44 ±5 

13.94 ±0.04 

13.88 ±0.05 

-55, ±40 

Fell 

1055.262 

0.898 

54 ±4 

49 ±5 

14.92 ±0.03 

14.88 ± 0.04 

-90, ±15 

Fell 

1063.176 

1.765 

165 ±5 

160 ±6 

14.73 ±0.03 

14.68 ± 0.04 

-90, ±35 

Fell 

1096.877 

1.545 

e 

132 ±4 

e 

14.76 ± 0.02 

-75, ±25 

Fell 

1125.448 

1.255 

93 ±4 

93 ±4 

14.83 ±0.02 

14.83 ± 0.02 

-75, ±35 

Fell 

1127.098 

0.483 

14 ± 3 

13 ±3 

14.69 ±0.09 

14.65 ± 0.09 

-55, ±5 

Fell 

1142.366 

0.681 

24 ±4 

27 ±3 

14.75 ±0.08 

14.79 ± 0.05 

-50, ±20 

Fell 

1143.226 

1.306 

96 ±5 

90 ±4 

14.78 ±0.03 

14.77 ±0.02 

-80, ±25 

Fell 

1144.938 

1.978 

231 ±6 

235 ±4 

14.68 ±0.06 

14.74 ± 0.04 

-95, ±30 

Felll 

1122.524 

1.786 

115 ±4 

115 ±4 

14.39 ±0.02 

14.42 ± 0.02 

-90, ±15 


^Measurements are given for FUSE detectors 1 and 2. Values refer to measurements with LiFl/LiF2 
for A > 1000 A and SiCl/SiC2 for A < 1000. All upper limits are Sa estimates and were derived 
assuming the transition has a breadth equal to the observed breadths for transitions of similar species 
(see §4). Rest wavelengths are from D. Morton 1999, private communication. 

'’Line contaminated by detector fixed-pattern noise. 

'’The N II result is based on measurements of SiClA and SiC2B data. Due to the overlap of the 
interstellar line by a broad stellar feature, we have used a model atmosphere kindly calculated by P. 
Chayer 2001, private communication, to normalize the spectrum before measuring this line (see text). 
The model is an excellent ht to the stellar lines in this region. 

‘'The measurements of P II in the two detectors are not in agreement, likely due to detector effects. 
We adopt the LiF2 measurement, but do not place too much conhdence in it (see text). 

"Line lies near the edge of the detector and therefore no accurate measurement is possible. 
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Table 5. Adopted Interstellar Column Densities 


Species 

log 

b [km s-’]’‘ 

Method^ 

HI 

19.97 ±0.03 


d 

H 2 

< 14.35 


AOD® 

C II 

> 14.55 


AOD 

C III 

> 13.85 


AOD 

NI 

> 15.35 


AOD 

Nil 

> 14.52 


AOD 

0 I 

> 16.35 


AOD 

0 VI 

14.49 ± 0.03 


AOD 

Si II 

> 15.20 


AOD 

P II 

13.46 ±0.10 

> 9 . 9 ’ 

AOD,COG’ 

P III 

< 13.4 


AOD 

P V 

< 12.7 


AOD 

S III 

14.47 ± 0.03 


AOD 

S IV 

< 13.7 


AOD 

S VI 

< 13.4 


AOD 

Arl 

14.02 ± 0.05 

9.9 ±0.8 

COG 

Fell 

14.80 ± 0.05 

22.9 ± 1.0 

COG 

Felll 

14.42 ± 0.05 

22.9 ± 2.0® 

AOD,GOG 


‘‘Adopted column density. All upper limits are Scr 
estimates. 

’‘Doppler parameter derived via curve-of-growth 
analysis. This quantity is not given for species for 
which the column densities were derived via the ap¬ 
parent optical depth. 

'‘AOD: Apparent optical depth method. Column 
densities are derived from direct integrations of the 
apparent column density profiles. COG: Curve-of- 
growth fit to measured equivalent widths of the lines 
for each species. 

‘’H I column density derived from STIS echelle ob¬ 
servations of Lyman-a from Howk et al. 2003. 

‘’The limits on H 2 are derived in §3. 

^The adopted P II column density is an average of 
the values derived using apparent column density and 
curve of growth analyses of the LiF2A data. The lat¬ 
ter assumes a fo-value equivalent to that derived for 
Ar I. Given the potential systematic uncertainties, 
the quoted uncertainties represent the linear combi¬ 
nation of the uncertainties in each method. 

®The Fe III column density has been derived using 
a curve of growth with the &-value (with twice the 
uncertainties) derived for Fe II. The AOD method 
gives consistent results. 






Table 6. 


Kinematic Properties of Select Species 


Species 

(wlsr)’" 
[km s“^] 

{o-obs)'° 
[km s“^] 

{c^oY 

[km s“^] 


0 VI 

-26.3 ± 1.3 

32.9 ± 1.1 

31.8± 1.1 

1.01 

Si II 

-20.7 ±0.5 

18.2 ±0.4 

16.1 ±0.4 

1.55 

P II 

-22 ±4 

14 ±3 

11±3 

0.58 

Sill 

-26.8 ±2.1 

17.0 ± 1.7 

14.7± 1.7 

0.79 

Ar I 

-9.7 ±0.6 

13.7±0.5 

10.7±0.5 

0.94 

Fe II 

-23.4 ± 1.1 

14.2 ±0.9 

11.4±0.9 

0.80 

Fe III 

-30.8 ±0.8 

15.3 ±0.6 

12.7±0.6 

0.95 


“^Average of the central velocities determined by fitting 
Gaussian models to the Na{v) profiles. The uncertainties 
are statistical only (they do not include the relative uncer¬ 
tainties associated with the FUSE wavelength solution or 
the absolute correction to the LSR frame discussed in the 
text). 

Average of the fitted Gaussian dispersions. 

‘^Gaussian dispersions corrected for the contribution of 
the instrumental line spread function, assumed to be 
<7lsf ~ 8.5 km s“^. 

‘^Average peak optical depth of the transitions used in 
the Gaussian fits. 



Table 7. WNM Gas-phase Abundances 


Species 

log(A/H )0 + 12'’ 

[A/H]b 

CII 

8.52 

> -1.94 

NI 

7.93 

> -0.55 

01 

8.74 

> -0.36 

Si II 

7.55 

> -0.32 

P II 

5.56 

-0.07 ±0.10 

Ar I 

6.48 

-0.43 ±0.06'= 

Fe II 

7.50 

-0.64 ±0.06 


“^Solar system abundances adopted from 
Grevesse & Sauval 1998 with the following 
exceptions. Argon: We adopt the Sofia & 
Jenkins 1998 value, which is the average of 
solar and local B-star photospheric abun¬ 
dances. Oxygen and nitrogen: we adopt 
the new photospheric determinations of Hol- 
weger 2001 . 

^’Normalized gas-phase abundance, 
relative to hydrogen: [A/H] = 

log[iV(X)/iV(HI)] - log(A/7J)o. This 
treatment neglects ionization effects as well 
as the effects of the large beam used for 
determining A^(H I). 

'’The apparent deficiency of Ar I is likely 
caused by the over-ionization of Ar I rela¬ 
tive to H I. See Sofia & Jenkins 1998 for a 
discussion of this effect. 
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Table 8. Census of Phases Toward von Zeipel 1128 


Phase 

Scale Height 

Ref.^ 

\ogN(H)'^ 

Fraction'^ 

Tracer'^ 

Warm Neutral 

0.4 kpc® 

1,2 

19.97 

53% - 86% 

HI 

Warm Ionized 

1.0 kpc 

3,4 

19.1 - 19.9 

12% - 45% 

SIII 

Warm-Hot Ionized 

2.3 kpc 

5 

^ 18-4 

> 1% - 2% 

0 VI 


'^Reference for scale height determinations: (1) Savage & Massa 1987; (2) 
Diplas & Savage 1994; (3) Haffner et al. 1999; (4) Savage, Edgar, & Diplas 
1990; (5) Savage et al. 2002. 

^’Estimated total hydrogen column density associated with each phase of the 
ISM (see text). 

‘^Fraction of total hydrogen column density associated with each phase of the 
ISM. A range of values is given to account for the various ionization corrections 
assumed. 

‘^Tracer used to derive the hydrogen columns associated with each phase of the 
ISM. 

®The H I scale height given here corresponds to the extended thick disk dis¬ 
tribution rather than the thin disk distribution (see, e.g., Diplas & Savage 1994; 
Dickey & Lockman 1990). 
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Fig. 1.— Ground-based spectra of ionized and neutral hydrogen emission towards vZ 1128. The top panel 
is a WHAM spectrum of Ha emission taken with a 1° beam at ~ 12 km s“^ resolution (Haffner et al. 2002). 
The spectrum shown here is an average of all pointings within 1?5 of vZ 1128 and therefore probes the Ha 
emission of a region of the sky with radius ~ 2°. The middle spectrum is from the Leiden-Dwingeloo H I 
21-cm Survey (Hartmann & Burton 1997). This survey has a 30' beam and 1.0 km s“^ resolution. The 
bottom spectrum is an NRAO 140-ft H 1 21-cm spectrum from Danly et al. (1992). These data sample a 
beam size of 21' at 1.0 km s“^ resolution. 
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Fig. 2.— The FC/5F spectrum of vZ 1128. Prominent interstellar lines are marked. For wavelengths longward 
of 1000 A, the unmarked lines are all either stellar or blends of stellar and interstellar features. Shortward 
of 1000 A we have only marked interstellar transitions of H I, O I, and N I. We have also noted the positions 
of stellar He II transitions in this wavelength range below the spectrum in the top panel. 
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Fig. 3.— Normalized profiles of several interstellar transitions versus LSR velocity, the majority taken from 
FUSE LiFlA and LiFlB spectra. These data have been binned by three detector pixels to 20.1 mA, or ~ 6 
km s“^, per data point. Stellar features and stellar/interstellar blends are marked with an asterisk. Other 
interstellar transitions in the velocity range are marked. The C III and N II profiles shown are from the 
SiClB and SiClA detector segments, respectively. Each of these profiles has significant contamination from 
stellar absorption. We have used a stellar model to remove most of the overlapping stellar absorption from 
the N II profile (see text). 
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Fig. 4.— The Green Bank 140-ft H I profile (top) and the O VI apparent column density profile {bottom two 
panels) observed toward vZ 1128. The O VI Na{v) profiles are shown with the best-fit Gaussian overlaid 
{middle) and three models of the Na{v) distribution assuming the O VI participates in Galactic rotation 
{bottom). The best-fit Gaussian has a dispersion cr = 31.8 ± 1.1 (after removal of instrumental effects). The 
rotation models all assume the O VI is distributed with a scale height of 2.3 kpc. The three models differ 
in their assumed cloud velocity dispersion (using cr = 10, 20, and 30 km s“^) and have been scaled to best 
match the red wing of the O VI profile. 
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Fig. 5.— Normalized gas-phase abundances of species found in the warm neutral medium along the vZ 1128 
sight line (points). No corrections have been made for ionization effects. The hatched area represents the 
range of values commonly found along Galactic halo sight lines for silicon and iron. For nitrogen and oxygen, 
we use the la confidence intervals of the distribution of measured values from Meyer et al. (1997) and Andre 
et al. (2003), respectively. For argon the values represent the range measured by Sofia & Jenkins (1998). The 
values given for phosphorous, which has few measurements for low density environments, represent the full 
range of recent measurements: the lower point from Welty et al. (1999) and the upper point from Howk et al. 
(1999). We note that Jenkins, Savage, & Spitzer (1986) derive an average value of [P/H] « —0.2 for WNM 
sight lines studied with Copernicus. We assume solar system abundances from Grevesse & Sauval (1998), 
except for argon (6.48; from Sofia & Jenkins 1998) and nitrogen and oxygen (7.93 and 8.74, respectively; 
from Holweger 2001). 
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Fig. 6.— Apparent column density profiles of several ionic species observed toward vZ 1128 that trace the 
warm neutral (Ar I, Si II, Fe II), the warm ionized (Fe III, S III), and the hot ionized (O VI) gas toward 
vZ 1128. The two singly-ionized species may also contain a contribution from the warm ionized medium. It 
should be noted that the profiles of Ar I and Si II likely contain unresolved saturated structure that causes 
the Na{v) distribution to be different from the true column density distribution. The analysis of the other 
transitions suggests relatively small saturation effects. Dashed lines in several panels represent unrelated 
interloping absorption lines. 
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Fig. 7.— Apparent column density profiles of Fe III compared with Fe II, S III, and O VI. In each case 
the Fe III profile has been scaled to match the other species. The WNM and WIM species trace each other 
quite well (within the limitations of the FUSE resolution), though there may be a slight excess of Fe III near 
vlsr ~ — 50 km s“^ compared with the others. The Fe III and O VI have nearly the same central velocity, 
though the latter, which traces hotter gas, is much broader. 
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Fig. 8.— Absorption line profiles of O VI A1031.926, C II A1036.337, C III A977.020, H I A937.803, and H I 
A930.748. The H I profiles have been normalized by stellar model atmosphere to account for the strong, 
broad stellar photospheric H I absorption in the negative velocity wings of the interstellar profile. Residual 
broad absorption in these profiles may be in part a result of mismatches between the data and the stellar 
model. The vertical lines show the approximate breadth of the C II profile at half intensity. No significant 
high-velocity (|ulsr| ^ 125 km s“^) absorption is seen in any of these profiles, though the negative velocity 
portion of the H I, C II, and C III profiles are strongly contaminated by stellar absorption. The excess flux 
at the bottom of the H I 937.803 A profile is likely residual terrestrial airglow emission. 




















